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Abstract 

We investigate the modification of antikaons in isospin symmetric and asymmet- 
ric nuclear matter. Using the leading s-wave couplings of the SU(3) chiral meson- 
^sD ! baryon Lagrangian we solve the coupled channel kaon-nucleon scattering equation 

selfconsistently. The in-medium antikaon propagator is calculated for different den- 
sities and different proton/neutron ratios. The spectral function of the antikaon 
is found to be broadened strongly, and its mass is shifted downward significantly. 
However, comparing the effective in-medium mass of the antikaon to the relevant 
charge chemical potential of neutron star matter, we find that kaon condensation is 
unlikely to occur in neutron stars. 
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Understanding the properties of strongly interacting matter under extreme conditions 
such as high temperatures or densities is of considerable interest in hadronic physics. In 
this context strange particles play a key role. In ultrarelativistic heavy-ion collisions, 
primarily exploring the high-temperature regime of the QCD phase diagram, strangeness 
enhancement is considered to be one of the possible signals for the formation of the quark- 
gluon plasma jTj. At low temperatures but high densities, as present, e.g., in the interiors 
of neutron stars, strangeness might show up in the form of a condensate of negative kaons. 
This possibility was first pointed out by Kaplan and Nelson E] and has received great 
attention since. The appearance of an antikaon condensate would soften the equation of 
state of the star, allowing more compact stellar objects. 

The basic considerations leading to this idea are relatively simple j3] : Since kaons, as 
the lightest mesons with strangeness, are the (pseudo-) Goldstone bosons of SU (3) chiral 
symmetry breaking, their interactions can be studied within chiral perturbation theory 
(xPT). At lowest order, the corresponding kaon-nucleon scattering Lagrangian contains 
two s-wave interaction terms: a constant scalar term due to explicit chiral symmetry 
breaking and a momentum-dependent vector term, the so-called Weinberg-Tomozawa 
term jHUHl- The scalar interaction is attractive for both, kaons and antikaons, while the 
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vector interaction is repulsive for kaons (K + and K°) but attractive for antikaons (K~ and 
K°). As a result, one expects an interaction which at tree level is weakly repulsive for kaons 
but strongly attractive for antikaons. Hence, despite large quantitative uncertainties, the 
antikaon mass should decrease in matter, eventually leading to kaon condensation, once 
the mass drops below the chemical potential of the surrounding leptons. 

Unfortunately, the situation is not that simple. Firstly, it should be noted that there 
is an ongoing controversial discussion about possible ambiguities related to the off-shell 
extrapolation from the .fTA-scattering amplitude to the kaon selfenergy [El HE FTP] . 
Moreover, the empirical value of the s-wave K~ nucleon scattering length is repulsive at 
threshold (JR,ea K ~ p = —0.78 fm). This discrepancy between experiment and the above 
considerations can be attributed to the existence of the A(1405) resonance just below the 
KN threshold, which gives rise to a repulsive contribution to the scattering amplitude 
at threshold. Apart from a few exceptions where the A(1405) has been introduced as 
an elementary field [TT], in most models it is dynamically generated through the KN - 
scattering process [121 EH UH HHj • In this case it cannot be treated perturbatively. 

On the other hand, there are various experimental indications that antikaons, once 
they propagate through dense nuclear matter, indeed feel an attractive potential. Studies 
of the energy shifts and the widths of the lowest levels in kaonic atoms probe the behavior 
of antikaons in nuclear matter of densities between zero and po- A collection of data 
over a wide range of atoms was analyzed by Friedman ^Sj. These results require an 
attractive potential for the antikaons, although different theoretical approaches lead to 
quite different values for the depth of this potential [TT| ITSj. 

This finding is supported by the kaon data in heavy-ion collisions. The azimuthal 
emission patterns of K + and K~ obtained by the FOPI and KaoS collaborations pHl EH] 
indicate a different in-medium behavior for kaons and antikaons, corresponding to repul- 
sive and attractive interactions, respectively. Naively, this also gives a simple explanation 
for the enhanced K~ / K + ratios which are found in nucleus-nucleus collisions as compared 
with (anti-) kaon production in NN: Lowering the K~ mass would reduce its production 
threshold and thus increase the yield [2U [221 12H] • However, the situation is complicated 
by the fact that in dense matter the K~ yield is enhanced, even without medium modi- 
fication of its mass. This is because there are hyperons around (e.g., owing to prior K + 
production) which enable secondary production mechanisms, as irY — > K~N . Hence, for 
a consistent interpretation of the data, we do not only need a good understanding of the 
in-medium modification of the (anti-) kaons, but also of the above production process. 

It is exactly this coupling of ixY to KN that turns out to be the key to the interaction 
of the K with matter, giving rise to the A(1405) resonance. This situation requires the 
treatment in terms of coupled channels, leading to the A(1405) resonance and connected 
effects. Performing the coupled channel calculation at finite density results in a A (1405) 
that is shifted upwards in energy [T2l ITSj [2U [23 • At lower energies around the original 
(vacuum) threshold the optical potential for the K~ becomes attractive, which would 
explain the data from kaonic atoms. 

Still, one has to be careful with this kind of analysis. In the first works that went 
beyond a mean field description to include the A(1405) (e.g. [121 EH 1211), the repercus- 
sion of the dynamics of the A(1405) on the antikaons were not included selfconsistently. 
However, this turned out to be of great importance: The selfconsistent A(1405) does not 
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move very much fTJEH^Hl, rendering the argument for the kaonic atoms questionable. 
Obviously, a comprehensive treatment of all the effects becomes necessary. This is what 
we aim at in this work. To this end, we basically carry out the following program, which 
is sketched in Fig. [U 
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Figure 1: Illustration of the selfconsistency scheme: Scattering equation for the antikaon- 
nucleon T matrix (upper line), in-medium antikaon self-energy (center), and dressed kaon 
propagator (lower line). 



In a first step we construct the KN T-matrix, taking into account the coupling of 
KN to all other important channels. Here the Lagrangian of chiral perturbation the- 
ory provides the basic interaction, i.e., the scattering kernel for the various mesons and 
baryons we take into account. As already indicated, the A (1405) appears slightly (about 
30 MeV) below the KN threshold which prohibits a purely perturbative treatment of the 
scattering process. Instead, we take the simplest diagram formed from C x pt as the kernel 
of the Bethe-Salpeter scattering equation, sum it to all orders and get the A(1405) as 
a dynamically generated resonance. In this context the coupling to other meson-baryon 
channels, in particular the 7rE channel, which opens about 70 MeV below the resonance, 
is crucial. 

Having constructed the scattering amplitudes in vacuum, we take the calculation to the 
dense medium, replacing the vacuum propagators by in-medium ones. The in-medium 
self-energy of the antikaon can be obtained from the T-matrix by closing a nucleon (hole) 
loop (see Fig. [TJ) . This leads to a modified kaon propagator which, in turn, should be used 
to calculate the T-matrix. This procedure has to be iterated until a selfconsistent result 
is achieved. 

A very important aspect in our calculations is the consideration of the momentum 
dependence of all calculated quantities. Calculating the antikaon selfenergy from the 
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in-medium scattering T-matrix means an integration over momentum, so the T-matrix, 
depending on energy and momentum separately in the medium, has to be known as a 
function of energy and momentum. This requirement becomes even stronger once we 
iterate the procedure to selfconsistency. 

This is not just a technical issue. In general the existence of the nuclear matter 
rest frame as a preferred frame of reference forces us to specify explicitly the energy 
and momentum dependence of in-medium properties. In addition, considering heavy- 
ion collisions, the optical potential of the antikaon is probed at finite momentum. At 
a moderate temperature of T = 80 MeV jTHj the antikaon has an average momentum 
of more than 300 MeV with respect to the matter rest frame. Momenta above ~ 300 
MeV are typical in the experimental data on heavy ion collisions. Nevertheless, in many 
calculations it has been assumed that the attractive potential for the antikaons, i.e. their 
change in the mass is independent of momentum f3J|2I|. This is not justified. To the 
contrary, a considerable momentum dependence of all in-medium properties is found |25j . 
see also section 

In this article we focus on s-wave interactions. Kaon condensation is thus connected to 
the propagator at zero momentum, p-wave antikaon condensation has been investigated 
elsewhere j2H I2H1 1201 EH EH EH) • One aspect here is the effect of hyperon-nucleon-hole 
contributions on the antikaon selfenergy. We briefly investigate these contributions in 
Sect. 13.41 As they turn out to be very small, but would lead to considerable technical 
complications, we carry out the selfconsistency problem with the s-wave part alone. 

The remainder of this work is organized as follows. In Sect. HI we briefly discuss 
the formalism involved in calculating the 577(3) meson-baryon scattering in a coupled 
channels scheme, based on the xPT Lagrangian. After that, in Sect.0 the model is used 
to investigate the behavior of antikaons in dense isospin symmetric matter. Among other 
aspects we study the effect of the selfconsistent treatment, the momentum and density 
dependence of the propagator, and the role of the in-medium modification of pions. In 
Sect.|3]the analysis is extended to isospin asymmetric matter, i.e., to systems with different 
densities of protons and neutrons, as found in realistic systems. The implications of these 
results for the question of antikaon condensation in neutron stars are discussed Sect. |5] 
Our main results are summarized in Sect. |UJ 

2 Meson— baryon scattering 

As outlined above (see Fig. [TJ) , the first step in our procedure to construct the in-medium 
kaon propagator consists in a coupled-channel calculation of the KN — > KN scattering 
amplitude, 



(K(k')N(p')\T\K(k)N(p)) = (2iT) 4 6 4 (k'+p'-k-p)u( P ')T KN ^ KN (k',p';k,p)u(p). (1) 



The T matrix elements are given by the Bethe-Salpeter equation [3*4"j . 
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which is graphically sketched in the upper line of Fig. ^ We will refer to this equation also 
as the T-matrix equation. The indices % and / denote the incoming and outgoing meson- 
baryon channels, while the index c indicates the intermediate states, i.e., the meson- 
baryon channels that form the loop. For brevity we have suppressed all momentum 
arguments and only kept the loop momentum I. 

Input to the T-matrix equation are the interaction kernel V, which will be specified 
in Sect. 12.11 and the two-particle propagator Q BS . The latter is basically the prod- 
uct of the single-particle propagators of the baryon B c and the meson M c in the loop, 
Qf s = iSb c Gm c - In this work we include all channels of the SU(3) meson and baryon 
octets coupled to total strangeness —1, dropping the S-channels that open only at larger 
energies. 



2.1 Interaction 

As discussed earlier, the interaction kernel V which enters the T-matrix equation is 
derived from the 577(3) chiral Lagrangian. The basic building blocks are the matrices 
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corresponding to the pseudoscalar meson and baryon octets, respectively. 

At lowest chiral order, the Lagrangian contains two meson-baryon interaction terms: 
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The first term, 
Cwt 



itr(B-f\r„B]) , with r A 



ud^u 



(5) 



(6) 



is the already mentioned Weinberg-Tomozawa term, which corresponds to s-wave inter- 
action. As pointed out in the Introduction, this term is crucial for a realistic description 
of KN scattering and will give a main contribution to the interaction kernel. 

The second term, C pw , corresponds to p-wave interactions. As we will show in 
Sect. 13.41 . it has only little effect in the KN sector. Of course, p-wave interactions must 
be included to describe medium modifications of pions. However, since pions are not 
selfconsistently coupled to the kaon sector in our model, they can be treated separately. 
This will be briefly discussed in Sect. 12.41 

At second chiral order there are three interaction terms due to explicit chiral symmetry 
breaking through the quark mass matrix M. = diag{m ui m^m^). These are the "Sigma 
terms" [SHI ESE EZ] , which are given by 



£x = b D tr(B{ X+ ,B}) + b F tr(B[ X+ ,B}) + b tr{BB)tr{ X , 



(7) 
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with 

X+ = 2B {u ] Mu ] + uMu) . (8) 

Although of higher chiral order, they are comparable in strength to the Weinberg-Tomozawa 
contribution when evaluated at the respective meson-baryon thresholds. Therefore they 
should be taken into account. 

While B can be related to the chiral condensate via the Gell-Mann-Oakes-Renner 
relation, B = —{qq)/f%, the constants bc,b F and b are examples of the so called low- 
energy-constants (LECs). These parameters are not constrained by symmetries and have 
to be fixed from phenomenology. One can use the predictions of the theory for the baryon 
masses to pin down the values of at least bo and bp- The values used for the calculations 
in this work are taken from |38| : 

b D = 0.061 GeV -1 , b F = 0.195 GeV -1 , b = -0.346 GeV" 1 . (9) 
2.2 Factorization of the T— matrix equation 

The difficult point in solving the Bethe-Salpeter equation, Eq. (J2J), is any momentum 
dependence in the vertices. For constant interaction kernels V, the equation can be 
factorized and reduces to a simple matrix equation which can straightforwardly be solved 
by inversion: 

T = V + VJT => T = (1-VJ)- 1 V. (10) 

Here J is a diagonal matrix in the space of meson-baryon channels c corresponding to 
the loop integrals 

f d 4 l 

J c = i J J^yS Bc G Mc ■ (11) 

It depends on the total 4-momentum of the pair only and, hence, as long as V is constant, 
the T matrix also depends only on the total 4-momentum of the meson-baryon system. 

However, in our case, V is not constant, because the Weinberg-Tomozawa term, 
Eq. (JOJ), leads to 4-momentum dependent vertices via the derivatives 9 M . To avoid this 
difficulty we follow the commonly used strategy which involves the following approxima- 
tions: First, we evaluate Cwt in the baryon rest frame pl%l ll5|. This amounts to replacing 

C WT -> itr(B^[F°,B}) , (12) 

which removes the 3-momentum dependence. Still, the vertices depend on the energies 
k° and k'° of the incoming and outgoing meson fields, respectively, i.e. 

V WT oc ±- (k° + k'°) . (13) 

Therefore we solve the T-matrix equation by applying the so-called on-shell factorization 
[3*%] . In practical terms this means that the T-matrix is calculated as in Eq. (jlOj) . and 
only in the final expression V is taken to depend on the energies k° and k'° of the external 
mesons. This approximation can be justified by the observation that any effect of the 
energy dependence of the vertices inside the loops can be absorbed in renormalization 
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constants and does not influence the on-shell part of the T-matrix which is related to the 
physical amplitudes, Eq. [TJ For a more detailed discussion of these issues, see Refs. [7J 

EiiiniiHHiiHaiini. 

As a result of the above procedure, the problem reduces to the multiplication and 
inversion of matrices spanned by the various meson-baryon channels. Conceptually, this 
is of course much easier than solving the original Bethe-Salpeter equation. Nevertheless, 
owing to the coupling of the various meson-baryon channels, the expressions for the T- 
matrix elements in terms of the vertex- and loop matrices V and J comprise several 
thousand terms each, and thus will not be given here explicitly. 

2.3 Regular izat ion 

We regularize the loop integral by employing twice subtracted dispersion relations in 
the calculation of its real part. The subtraction parameters can be fixed by fitting the 
resulting scattering amplitudes to the empirical scattering lengths [4*T| I42j . 

This approach has some advantage over the use of cut-offs or form factors: On the one 
hand, it keeps the vacuum-amplitude Lorentz-invariant. On the other hand it helps to 
overcome a technical problem of the selfconsistency procedure we aim at: Iteratively using 
the scattering matrix and the modified antikaon propagator as mutual input (see Fig. HJ 
requires the knowledge of these functions on infinite energy and momentum intervals. At 
large momenta, the antikaon propagator can be approximated by the free one. However, 
it is not immediately clear how to treat the scattering amplitude in that respect. Use of 
a cut-off will introduce artificial peaks into the T-matrix, especially at higher energies 
and momenta. This is not the case for the amplitude obtained by a twice subtracted 
dispersion relation, which is smooth enough to be restricted to a finite interval as well. 

We find that our vacuum scattering amplitudes are in good agreement with the KN 
amplitudes of [5%] . 

2.4 Pion— nucleon interaction 

The formulation of the KN scattering problem in terms of coupled meson-baryon chan- 
nels is the key ingredient of our model. In particular, the A(1405) will only ciriSG clS cL 
dynamically generated resonance if the KN channel is coupled to the 7rE channel. Con- 
sequently, the in-medium amplitude will feel any medium-modifications of the pion. 

The behavior of the pion in dense matter has been studied in great detail. The main 
sources of the pion selfenergy are particle-hole and delta (A(1232))-hole excitations. To 
capture these and other important features, in the pion sector we do not stick to the 
Lagrangian set up in Sect. 12.11 but include additional phenomenological terms which are 
taken from the literature. Here we basically follow Ref. |4"5] . 

The p-wave pion-nucleon interaction can be obtained from C pw in Eq. (jSJ. It is given 

by 

AriV = ^ 7 Vr^ • dj, (14) 
m w 

with the pseudovector coupling constant /j\r = 1.01. 
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The coupling to the A looks quite similar: 



C vNA = -—^^a-d^f + h. C. 



(15) 



Both, the irNN and 7riVA vertices, are supplemented by a monopole form factor 



with A = 550MeV. Finally, short-range correlations are taken into account via Migdal 
parameters g' NN = 0.8 and g' NA = g' AA = 0.5. No indications for p-wave pion condensation 
are observed for the densities considered here (up to 5po)- 

Having calculated the pion propagator in dense matter, this propagator is put into 
the loops with the A and E baryons. These loops are then coupled to the other meson- 
baryon channels to obtain the T-matrix element. As pointed out before, we do not further 
modify the pion itself, i.e., we do not aim for selfconsistency in the case of the pion. This 
is justified because the above model of the pion [131 EH ESI EE] is already designed to give 
a realistic result as it stands, without further iteration. 

3 The kaon propagator in symmetric nuclear matter 

We can now proceed to calculate the in-medium antikaon propagator following the self- 
consistent scheme outlined in the Fig. ^ 

A typical example for the resulting antikaon selfenergy is given in the upper panel 
of Fig. |2l where the real and imaginary parts of at normal nuclear matter density, 
P — Po — 0.16 fm -3 , are displayed as functions of energy at a finite three-momentum of 
100 MeV. The corresponding antikaon propagator Gk is shown in the lower panel. The 
figure displays the typical properties of the antikaon when it is modified by the scattering 
processes in the medium. 

The main peak of the spectral function (the imaginary part of Gk) is shifted to lower 
energies compared with its vacuum position at uj = a/ m 2 K + k 2 ~ 505 MeV. However, 
this mass shift is not dramatic. A much stronger effect is the broadening. Taking the full 
width at half maximum, the antikaon is "spread" over roughly 130 MeV. 

In this section we want to investigate in some more details which are the main physical 
sources for this behavior and how this depends on the kaon-momentum and on the den- 
sity of the surrounding medium. We thereby concentrate on symmetric nuclear matter. 
Asymmetric matter will be discussed in Sect. El 

3.1 Effect of selfconsistency 

In this subsection we briefly discuss the relevance of finding a selfconsistent solution, rather 
than dressing the kaon perturbatively. Technically, once the selfconsistency calculation is 
set up and stable, the procedure works very well, converging after 3-4 iterations. 

The effect of the selfconsistency iterations on the KN scattering amplitude is shown 
in Fig. El where the imaginary parts of Tkn for vanishing three-momentum and densities 
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Figure 2: Real part (solid) and imaginary part (dotted) of the antikaon selfenergy (upper 
panel) and antikaon propagator (lower panel) as functions of energy w at a three-momentum of 
\k\ = 100 MeV and a density of p = pq. 
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p = Po (upper panel) and 5po (lower panel) are plotted as a functions of the total energy. 
For comparison we also show the vacuum result (dotted lines). Here one can clearly see 
the KN threshold at 1.43 GeV and the A(1405) resonance below threshold. When the 
density effects are included by taking into account the Pauli blocking of the nucleon but 
using the unmodified kaon propagator one arrives at the dashed lines. Because of the Pauli 
blocking the KN threshold is shifted to higher energies and, as a result, the A(1405) is 
also pushed upwards. 

However, in the selfconsistent calculation (solid lines in Fig. |3J), the change in the 
antikaon mass and especially the broadening of the kaon spectral function (see below) 
re-enables scattering into lower lying states and the resonance is pulled down again. This 
behavior has been discussed first in Ref. [HI 1213 • 
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Figure 3: Imaginary part of the T-matrix QmTj^ N at = as a function of energy. Upper 
panel: p = po, lower panel: p = 5po- Dotted lines: vacuum result, dashed lines: in-medium 
result, only Pauli-blocking, solid lines: result after selfconsistency is reached. 
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Fig. |3] shows the corresponding effect of the selfconsistency iteration on the antikaon 
propagator (imaginary part) at p = p and two different 3-momenta. Again, the dashed 
lines, labeled "Pauli-blocked", mark the results after the first iteration step, i.e., when only 
the Pauli-blocked T-matrix elements (cf . dashed lines in Fig. EJ) are used in the selfenergy 
integrals. The selfconsistency iterations broaden the antikaon and shift it further down 
(solid line). 

Also note that in the "Pauli-blocked" calculation there is a double peak structure 
in QmG-K which becomes more prominent at higher momentum. This feature is largely 
washed out by the selfconsistency iterations. 
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Figure 4: Imaginary part of the antikaon propagator ^rnGx at normal nuclear density p = po 
as functions of energy. Upper panel: \k\ = 100 MeV, lower panel: \k\ = 300 MeV. Dashed lines: 
after one iteration, solid lines: after selfconsistency is reached. 
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3.2 Momentum and density dependence of the antikaon propa- 
gator 

Fig. |3] already indicates the structure of the antikaon spectral function as a function of 
momentum. This is summarized in Fig. |5] which shows the imaginary part of the kaon 
propagator at density p = po for three different 3-momenta, \k\ = 100 MeV, 200 MeV, 
and 300 MeV. 

As mentioned above, one typically finds two peaks. It turns out that for low values of 

— * — * 

| A; | most of the strength is contained in the lower peak (solid line). With increasing \k\, 
this peak becomes gradually reduced while the upper peak grows (dashed line). Finally, 
at large momenta, most of the strength is found in the upper peak (dotted line). 
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Figure 5: Imaginary part of the selfconsistent antikaon propagator at density p = po as function 
of energy. The various curves correspond to three different 3-momenta: \k\ = 100 MeV (solid), 
200 MeV (dashed), and 300 MeV (dotted). 

The nature of the two-peak structure can be easily understood in terms of a simple 
two-level model, describing the coupling of a bare antikaon state to a A(1405)-hole state. 
At low momenta the unperturbed energies of these states are very close to each other, the 
A(1405)-hole state being only about 30 MeV below the bare kaon. However, turning on 
the interaction, the splitting of the peaks becomes larger through level repulsion, which 
shifts the "antikaon peak" somewhat up and the "A(1405)-hole peak" further down in 
energy. The latter receives most of the strength at low momenta while at large momenta 
the kaon branch dominates. This is quite analogous to the well-known case of a pion 
coupled to a A-hole state |43j . However, it should be noted that this picture is rather 
crude since, as a consequence of our coupled-channel approach, we cannot resolve the 
peaks in the propagator into particular quasiparticle states. 

The above features of the kaon propagator are subject to a strong variation with 
density. This is illustrated in Fig. |H1 where the imaginary part of the kaon propagator 
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is displayed as a function of energy at fixed 3-momentum \k\ = 300 MeV and the three 
densities p = p , 2 p , and 5 p . Obviously, the strength that at p = p is found in the 
higher-lying peak is still in the lower branch at the higher densities. This behavior is 
consistent: At higher densities the medium modifications are stronger. Thus the effect of 
increasing 3-momentum as discussed above sets in only at even higher momenta. 
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Figure 6: Imaginary part of the selfconsistent antikaon propagator at 3-momentum \k\ = 
300 MeV. The various curves correspond to three different densities: p = po (solid), 2po (dashed), 
and 5/?o (dotted). 

At very high momenta, the free dispersion relation of the antikaon is regained. This 
provides some technical help, allowing the use of the free antikaon beyond oj — 1.5 GeV 
and \k\ = 1.2 GeV. 

3.3 In— medium modification of the pions 

The importance of the pion in the KN scattering processes was already emphasized. As 
pointed out earlier, the pion undergoes considerable modifications when placed in dense 
matter. In the following, we want to study the effect of these medium modifications of 
the pion propagator on the kaon propagator through the coupled-channel mechanism. 

Our results are summarized in Fig. The pion propagator at normal nuclear matter 
density is displayed in the upper panel. At the indicated pion 3-momentum of 300 MeV, 
the broadening effect is maximal. 

Typically, the broadening of the pion in dense matter has some bearing on the 7T- 
baryon loop function in the form of a bump around threshold. As an example, the central 
panel of Fig. shows the imaginary part of the 7rS channel - loop function, 3m4s, at 
normal nuclear matter density and three-momentum \q\ = 100 MeV. The effect seems to 
be sizeable, the difference to the vacuum curve is quite pronounced. 
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Figure 7: Upper panel: Real part (dotted line) and imaginary part (solid line) of the pion 
propagator at | k\ = 300 MeV as functions of energy. 

Central panel: Imaginary part of the irY< loop function at \q\ = 100 MeV as function of the 
invariant mass, incorporating dressed (dashed line) and bare (solid line) pion propagators. 
Lower panel: Imaginary part of the antikaon propagator at \k\ = 100 MeV as function of energy, 
incorporating dressed (solid line) and bare (dashed line) pions. 
All calculations in this figure have been performed for p = pq. 
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On the other hand, the lower panel indicates the difference it makes for the imaginary 
part of the kaon propagator, QmG-^, whether the bare or the dressed pion propagator is 
used as input. The figure shows the self consistent result, again at normal nuclear matter 
density and \k\ = 100 MeV. Obviously, very little is left of the effect of pion dressing, 
in agreement with while in ^Hj a stronger effect has been found. The reason is 
the strong momentum dependence of the pion selfenergy. At the \k\ = 300 MeV as 
in the upper panel, its effect is most pronounced. However, already the loop function 
(central panel) is obtained by integrating the in-medium pion over all momenta. The 
kaon propagator (lower panel) is the result of yet another integration - not counting the 
various iterations to achieve a selfconsistent result. 

3.4 p— wave antikaon selfenergy 

As mentioned earlier, we neglect contributions to the antikaon selfenergy arising from 
the p-wave vertices of Eq. (JHJ) in our model. The reason is that these contributions are 
expected to be small, but would considerably complicate the selfconsistency problem. 

To check the (un-) importance of the p-wave vertices for the kaon propagator, we 
include their contributions to the selfenergy, but only in a perturbative way. This means, 
we first calculate the selfconsistent kaon propagator obtained with s-wave interactions 
only, as before. Then we add the p-wave contributions to the selfenergy, i.e., A-hole and 
X-hole diagrams as depicted in Fig. |S1 As in the pion case we use a monopole form factor 
with A = 550 MeV. For the Migdal parameters we take the classical value of g' = 1/3. 
The resulting total selfenergy is then used to calculate the modified kaon propagator, but 
the procedure is not iterated again. 




Figure 8: Diagrams corresponding to the p-wave contribution to the antikaon selfenergy. 



— * 

The result for p = p and \ k\ = 400 MeV is shown in Fig.El Although the chosen three- 
momentum roughly corresponds to the value where the p-wave effect has its maximum, 
it is obviously very small. This is due to the position of the p-wave selfenergy at very 
low energies. In fact, it lies entirely in the space-like region. Thus the kaon with its bare 
pole at uj ~ 600 MeV does not feel much of the p-wave selfenergy. Hence, even though 
somewhat larger effects have been found when the p-wave interactions are included in the 
meson-baryon scattering equation ^Hl 1213 EH] , their negligence seems to be justified. 
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and p-wave (solid) contributions. 



4 Asymmetric nuclear matter 

In this section we extend our discussion to isospin asymmetric matter, i.e., to the case 
of unequal densities of protons and neutrons. This case is relevant for the description of 
neutron star matter which is expected to consist of roughly 90% neutrons and 10% protons 
(see Sect. EJ). The calculation of the KN scattering and the antikaon selfenergy in the 
case of asymmetric nuclear matter follows the scheme developed in the previous sections. 
To discuss the systematic behavior we have investigated various proton-to-neutron ratios 
at various densities. In this context a somewhat more precise terminology is required: 
The density p referred to is always the total baryon number density. Since we consider 
the zero-temperature case, T = 0, the only baryons occurring in sizeable number are 
the nucleons. 1 Of these, a fraction x p consists of protons. The symmetric case is thus 
characterized by x p = 0.5. 

4.1 Kaons in asymmetric matter 

In Fig. EH the kaon spectral function at vanishing 3-momentum is shown for a proton frac- 
tion of 10% and three different densities. The most important difference to the symmetric 
case is the fact that, due to the isospin breaking surrounding, the propagators of K~ and 
K° are no longer identical, but the K° is much more affected by the dense medium. This 
can be explained by the fact that the Weinberg-Tomozawa term is strongest in the K°n 
and K~p channels. Therefore, since the neutron density is ten times higher than the pro- 
ton density, the K° undergoes a stronger modification, although the various contributions 

1 For a nudeon gas with 90% neutrons, the neutron Fermi energy reaches the mass of the (undressed) A 
only at p = 6.7po 
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get thoroughly mixed in the selfconsistency iterations. 

The variation of QruGx with the proton-neutron asymmetry at different densities can 
be studied by comparing the curves in Fig. ^2 This figure shows a smooth variation of 
the properties of the antikaons with density and isospin composition. 

Finally, Fig. El shows the imaginary part of the antikaon propagator, now comparing 
directly the different densities, at an asymmetry of x p = 0.1. As mentioned, the growing 
density affects the K° much more than the K~ . However, the overall effect is not dramatic, 
the peaks of the spectral functions being shifted by no more than ~ 200 MeV from the 
vacuum positions at 500 MeV. It also becomes obvious that the width of the antikaon 
shrinks again at high density. This is of course an effect of shrinking phasespace. 

4.2 Effect of s— wave pion— nucleon interactions 

In Sect. IH.HI we have seen that the medium modifications of pions have a rather small 
impact on the kaon propagator in symmetric nuclear matter. We have thereby restricted 
ourselves to include p-wave pion-nucleon interactions, as outlined in Sect. 12.41 However, 
there are also known s-wave interactions which we have neglected so far. The contribution 
of the Weinberg-Tomozawa term to the s-wave pion selfenergy is usually held respon- 
sible for the existence of deeply bound pionic states in heavy nuclei (1H1- The complete 
expressions for asymmetric matter in two-loop chiral perturbation theory are given by a 
calculation by Kaiser and Weise Their analytic expressions depend mainly on the 
difference of neutron and proton Fermi momenta, which means that the s-wave selfenergy 
is especially important in isospin asymmetric matter. 

It is therefore appropriate to check the influence of s-wave pion-nucleon interactions 
in the context of asymmetric nuclear matter. To that end we employ the results of 
Ref. jlHj to calculate the pion selfenergy. The effect on the 7r~-propagator is illustrated 
in the upper panels of Fig. EH for asymmetric nuclear matter with 10 % protons at a total 
density of 2p . We see that the n~ is pushed upwards in energy, i.e., the s-wave selfenergy 
is repulsive here. 

However, we find again that the impact of the s-wave pion selfenergy on the antikaon 
propagator is rather weak. The lower panels of Fig. 1131 show the K~ propagator after the 
first iteration when the pions in the underlying ttY loops are modified by p -wave and 
s + p -wave interactions. We infer from the figure that the s-wave 7r-selfenergy has some 
bearing on the antikaon propagator, but it is quite small. The zero of the real part of 
Gk- (the in-medium antikaon mass) is shifted upwards by just a few MeV (right panel 
in Fig.HHJ). 

Thus the pion modifications in general have only minimal impact on the in-medium 
kaon. 

5 The kaon mass under neutron star conditions 

As mentioned in the Introduction, particular interest in the properties of antikaons in 
dense nuclear matter is caused by the possibility of kaon condensation in neutron stars 

□ HO 
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Figure 12: Imaginary parts of the antikaon propagators as functions of energy at x p = 0.1, \k\ 
for different densities. K° on the left, K~ on the right. 



At not too high densities, neutron star matter consists of neutrons, protons, and 
electrons which are maintained in beta equilibrium through processes like n — > p + e~ + v e 
and p + e~ — > n + v e . For neutron stars older than a few minutes, the temperature is 
practically zero on the nuclear scale, and the neutrino mean free path is larger than the 
radius of the star, such that neutrinos can freely leave the system. Hence, the electron 
chemical potential is equal to the difference between the chemical potentials of neutrons 
and protons, p e = \i n —\i v . Moreover, neutron star matter has to be electrically neutral. In 
the above case, this means that the densities of electrons and protons must be equal, p e = 
p p = x p p. Employing realistic equations of state, one typically finds x p ~ 0.1 (see below). 
For the electron chemical potential, this roughly translates into ji e ~ 150 MeV (p/po) 1 ^ 3 - 

With increasing density, various other weak processes may become energetically pos- 
sible, leading to the occurrence of new particles. In fact, according to the estimate above, 
already at p = po, p e is well above the muon mass, and the conversion e~ — » pT +17^ + v e 
should lead to a finite density of muons. At higher densities, strangeness-changing re- 
actions should also occur. For instance, nucleons from the top of their Fermi seas could 
convert into hyperons. 

Similarly, electrons and muons could convert into kaons, 

— > K~ + v e , pT — > K~ + v e , (17) 

once the electron (= muon) chemical potential exceeds the in-medium mass of the K~ 
|501 IBTl 152*] . This is the case we are interested in. Obviously, the conversion of electrons 
and muons to kaons would reduce the number of fermions and thus their contribution to 
the pressure. On the other hand, the kaons might form a condensate at zero momentum 
that just provides a background of negative charge but, as a condensate of bosons, does 
not exert a degeneracy pressure. In addition, processes like n —>■ p+K~ would increase the 
proton fraction in the system. Since the isospin-asymmetry part of the nuclear interaction 
is repulsive, this would cause a further decrease of the pressure. Thus the formation of a 
condensate could considerably soften the nuclear equation of state "2HHS01- in Ref. joT"] . 
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Figure 13: Upper panel: Imaginary (left) and real part (right) of the tt~ propagator as functions 
of energy for \k\ = 300 MeV. Solid lines: results including s-wave and p-wave contributions, 
dashed lines: p-wsve only. 

Lower panel: Imaginary (left) and real part (right) of the K~ propagator after the first iteration 
as functions of energy at \k\ = 0, containing pions with including s-wave andp-wave selfenergies 
(solid line) and p-wave selfenergy only (dashed line). 

All calculations in this figure correspond to total baryon density p = 2po and proton fraction 
x p = 0.1. 
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this effect has been employed to predict a maximum mass for neutron stars (better: 
nucleon stars) of about 1.5 solar masses and, as a consequence, the existence of a large 
number of low-mass black holes in our galaxy. 

As pointed out above, a precondition for the conversion processes Eq. (j!7| ). is that 
the electron chemical potential must exceed the in-medium kaon mass, rn* K ^. Strictly 
speaking, the relevant quantity is not the mass, but the energy of the kaon. However, 
as we have seen in Sect. 13.41 p-wave contributions to the kaon self-energy are small, 
and therefore the lowest energies are found at vanishing 3-momenta where they can be 
identified with the mass. Since the vacuum mass is much too large to be reached by 
/i e at any realistic density, the question whether or not kaon condensation takes place in 
compact stars obviously depends crucially on how fast m* K „ drops as a function of density. 
In the following we want to investigate this point with our model. 

The mass of a given (quasi-) particle can be defined by the poles of the propagator. 
As explained above, we can restrict ourselves to the case of vanishing 3-momentum. We 
also ignore the imaginary part of the selfenergy at this stage and postpone its discussion 
to the end of this section. Then, the in-medium kaon mass is given by 



Alternatively, we could define the mass as the maximum of the spectral function. It turns 
out that the difference between these two definitions is small for the cases of interest. 

As we have seen in Sect. |3J the kaon propagator depends on both, total density p 
and proton fraction x p of the surrounding medium. Hence, in order to proceed, we have 
to know Xp clS £L function of p. For neutron star matter this is in principle fixed by 
the requirement of beta equilibrium and neutrality. However, to apply these conditions, 
we need to know the equation of state which determines the relation between chemical 
potentials and densities of the involved particles. As for the leptons, we can safely employ 
ideal gas relations. On the other hand, this does not work for the hadrons where an 
ideal gas equation of state yields unrealistically small proton fractions. Ultimately, of 
course, all hadrons should be described within a unified approach, where the in-medium 
pions and kaons are used to calculate the dressed nucleons and vice versa. Lacking such 
a description at the present stage, we have consulted the literature to get some insight: 
In Ref. [53], Akmal and collaborators compare x p in neutral beta-stable matter resulting 
from a variety of nucleon-nucleon interactions, such as the Bonn, Urbana and Argonne 
models as well as the Nijmegen results. Typically, the values are of the order x p = 0.1, with 
a general tendency that x p is rising with density. (For instance, taking their calculation 
based on the Argonne NN interaction vi$ and including boost corrections (£„) and a 
parameterization of three nucleon forces (UIX*), one finds x p = 0.06 at p = po and 
x p = 0.14 at p = 5p , while the Vi 8 interaction alone gives just x p = 0.095 at p = 5p .) 
In view of these and other uncertainties and the fact that our results are not extremely 
sensitive to small variations of x p , we take a fixed value of x p = 0.1 in our analysis. 

In Fig. we show the real part of the inverse kaon propagator at vanishing 
momentum 3-momentum as a function of energy. The three panels correspond to three 
different densities, p = p , p = 2p , and p = 5p - For completeness both the K~ (solid 
lines) and the K (dashed lines) are shown, although only the K~ mass is relevant for the 
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question of kaon condensation. For comparison, we also indicate the dispersion relation 
of the free kaon (dotted lines). 

The resulting effective kaon masses are listed in Table H As we have seen before, the 
neutral antikaons (which are irrelevant in the context of kaon condensation in neutron 
stars) are more strongly affected by medium effects than the negative ones. 

For comparison with the literature we also list the effective mass for symmetric mat- 
ter. At p = po we find m lf ymm = 438 MeV, corresponding to a optical potential depth of 

V Q p t = Ej^ N /2m^ ~ —60 MeV. This is in good agreement with [T7], where Vj£ t « —55 MeV 
has been extracted from an analysis of kaonic atoms. 

The masses are also plotted in Fig. ED as functions of p. The K~ masses are indicated 
by ▲ symbols. We find that the numerical results are well described by a linear fit 

m* K .{p) = m K -(p = 0)(l " c ^) (19) 

with c = 0.06 (solid line). 

To decide whether the mass drop is sufficient to trigger kaon condensation we need 
to know the corresponding values of the electron chemical potential /x e . To that end, 
we consider a gas of non-interacting electrons and muons and determine pi e by requiring 
charge neutrality together with the protons in the system. The result is indicated by the 
dashed line in Fig ED The corresponding numbers are also listed in Tabled 

Obviously, rn* K - is well above /i e , even at five times nuclear matter density. Even 
though the knowledge of nuclear interactions and the hadronic equation of state under 
neutron star conditions is not very good, the variations in x p and /i e caused by the 
uncertainties in present day calculations are not large enough to support the condensation 
scenario below 5po- Extrapolating our results to higher densities, we predict the onset 
of K~ condensation at about 8po- However, it is very unlikely that our model can be 
trusted up to such high densities. In the hadronic sector there are several competing 
channels which probably open earlier and which we have neglected so far. These could 
be, for instance, negative hyperons or pions. Note that each additional negative source 
decelerates the rise of p e and therefore shifts the crossing point with the kaon mass to 
higher densities. Eventually, a phase transition to quark matter should take place. Hence, 
although a more detailed investigation of all these possibilities is certainly necessary, we 
conclude that the occurrence of kaon condensation in neutron stars is rather unlikely. 



p/po 


p e [MeV] 


[MeV] 


m^o [MeV] 


rrq vmm [MeV] 


1 


141 


459 


420 


438 


2 


170 


430 


385 


407 


5 


221 


347 


294 


324 



Table 1: Total baryon density and corresponding values of the electron chemical potential and 
the in-medium masses of K~ and K° for x p = 0.1. Rightmost column: in-medium mass for 
x p = 0.5 
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Figure 14: Real part of the inverse kaon propagator at vanishing momentum 3-momentum as 
a function of energy for x p = 0.1 and three different densities: p = po (upper panel), p = 2po 
(center), and p = 5po (lower panel). Solid line: K~ , dotted line: K , dashed line: free kaon 
dispersion law, dashed-dotted line: K for x p = 0.5 
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Figure 15: In-medium K~ mass rn* K _ (A) and electron chemical potential jjL e (dashed line) 
versus baryon density at x p = 0.1. The solid line represents a linear fit through the calculated 
points. Also indicated are the K° mass (t) and the in-medium mass for the symmetric case 
(•)• 



Finally, we should recall that in our determination of the in-medium kaon "mass" 
(cf. Eq. (|18|l). we have neglected the imaginary part of the selfenergy. Therefore, one 
might raise the question whether the consideration of the kaon width could change our 
conclusions. However, this seems not to be the case: As one can see in Fig. ^1 even 
the low-energy tails of the spectral functions are far away from the values of /i e given in 
Tabled 

In this context we should emphasize that our present approach is not suited to describe 
the phase transition to a kaon condensed phase itself. If we further increase the density, 
eventually part of the spectral function will move below /i e while another part still lies 
above, and it would be unclear whether or not the conversion process Eq. (117)1 takes 
place 2 . This is an artifact of our model where we calculate the modified kaon propagator, 
leaving the ground state of the system, i.e., the Fermi sea of nucleons, unchanged. It is 
clear that, once kaon condensation sets in, the kaons with vanishing 3-momentum are 
part of the ground state, and the width of the propagator has to vanish at u = fi e . Hence, 
our present approach is only justified as long as the modifications of the nucleonic ground 
state are small in comparison with the kaon sector, i.e., as long as we are far away from 
the phase transition. In this sense, at p = 5po we are still on the safe side. 

2 Strictly speaking, this problem exists already at lower densities in our model: If we cut the kaon 
selfenergy in Fig. ^ we see that there is a contribution corresponding to the decay K~ — ► 7r°Ap _1 , 
where p~ x denotes a hole state in the Fermi sea of the protons. Moreover, the pion can further decay 
into a nucleon-hole state. Therefore the imaginary part of the kaon propagator opens at an energy 
approximately equal to the A-proton mass difference, where we have neglected the Fermi energy of the 
protons. This means that part of the kaon spectral function is formally below the electron chemical 
potential as soon as the latter exceeds this value which is the case at about 2p . 
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6 Summary 



We have investigated the properties of antikaons in dense nuclear matter, focusing on the 
regime of zero temperature and high densities. Taking the most important KN vertices 
derived in xPT: the Weinberg-Tomozawa term and the a terms, we have constructed the 
T-matrix for meson-baryon scattering within a coupled-channel approach, thereby dy- 
namically generating the A(1405) resonance. Regularization of the scattering amplitudes 
was accomplished by means of twice subtracted dispersion relations. The values of the 
subtraction constants have been fixed by a fit to the experimental scattering lengths. 

Having developed a satisfactory description of the scattering process in vacuum, the 
calculations have been extended to the dense medium. The medium was treated in the 
nucleon gas approximation, i.e., assuming a free Fermi gas of protons and neutrons. This 
leads to a modified T-matrix through Pauli blocking the occupied nucleon states. On the 
other hand, closing a nucleon (hole) loop in the T-matrix, we obtained the in-medium 
self-energy of the antikaon. In turn, the resulting kaon propagator was used to (re-) calcu- 
late the T-matrix. This procedure was iterated until a selfconsistent result was achieved. 
Convergence was reached after some four to five iterations. The actual implementation 
of the selfconsistency scheme requires the determination of the scattering amplitudes and 
the antikaon propagator for the entire energy-momentum plane. A suitable method to 
implement this requirement was developed. 

We first applied the model to isospin symmetric nuclear matter. We found that the 
A(1405) resonance is moved to higher energies and becomes strongly broadened. The 
antikaon also receives a finite width and is mostly shifted to lower energies. Its spec- 
tral function is found to be strongly momentum dependent, emphasizing the need for a 
calculation that correctly incorporates the finite 3-momentum of the scattering process 
and the in-medium antikaon. The spectral function at zero momentum is broad but the 
antikaon mass still seems to be well defined. At finite but small momenta there are two 
branches that stem from the in-medium excitations at lower energies and the antikaon 
pole at « 500 MeV. At higher momenta the antikaon dispersion relation approaches 
its free form. 

Besides the momentum dependence, we found the selfconsistent treatment of T-matrix 
and kaon propagator to have the largest effect on the results, while other features, like 
the in-medium modification of the pions and p-wave K AMnteractions were found to be 
less important. 

We then applied our model to isospin-asymmetric matter as found inside neutron 
stars. Although the general dependence on density and momentum is similar to that seen 
for symmetric matter, the most important difference is that, in asymmetric matter, K~ 
and K° develop distinct spectral functions. This could be traced back to the Weinberg- 
Tomozawa term which is strongest in the K°n and K~p channels. 

Finally, we investigated the possibility of K~ condensation in neutron stars. Conden- 
sation should set in when the K~ mass falls below the electron chemical potential in elec- 
trically neutral beta equilibrated matter. Neglecting the width of the K~ , we found that 
the mass can be described rather well by a linear fit, m* K _(p) = mx~{p = 0)(1 — 0.06p/po)- 
This drop is too slow to support kaon condensation for densities up to at least 5po- 
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